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ABSTRACT 


A  study  of  possible  configurations  for  a  Large  Aperture  Seismic  Array  has  been 
completed.  An  array  having  21  identical  subarrays  located  on  three  concentric  circles 
has  been  found  to  yield  the  most  satisfactory  pattern  in  wave  number  space  of  all  the 
configurations  tested.  Patterns  for  some  alternative  placements  of  subarrays,  including 
that  of  the  experimental  LASA  in  Montana,  have  been  included  in  this  report. 

This  study  of  patterns  in  wave  number  space  has  vieided  the  suggestion  that  a 
LASA  having  a  diameter  of  200  km  should  be  composed  of  subarrays  from  10  to  15  km 
in  diameter.  Such  an  increase  of  subarray  size  above  the  7  km  diameter  subarrays  of 
the  experimental  LASA  would  require  the  use  of  less  regular  subarray  geometries  than 
those  which  have  been  used  in  Montana. 

A  sensitivity  function  for  patterns  has  been  developed.  This  function  can  be  used 
to  predict  the  change  in  patterns  which  might  result  from  changes  in  seismometer  or 
subarray  positions.  Since  the  sensitivity  function  predicts  possible  changes  in  patterns, 
it  can  be  used  to  set  bounds  upon  changes  in  an  array  which  can  be  made  without  severe¬ 
ly  changing  the  pattern.  The  tight  bounds  imposed  by  the  sensitivity  function  can  be  re¬ 
laxed  if  the  pattern  resulting  from  any  anticipated  change  in  position  is  actually  computed. 
Since  very  little  computer  time  is  required  to  compute  a  pattern,  this  procedure  is 
highly  recommended. 


Accepted  for  the  Air  Force 
Stanley  J.  Wisniewski 
Lt  Colonel,  USAF 
Chief,  Lincoln  Laboratory  Office 
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"We  will  have  rings  and  things  and  fine  array." 
Shakespeare,  Taming  of  the  Shrew,  Act  II, 

Scene  1. 


I.  INTRODUCTION 

An  experimental  Large  Aperture  Seismic  Array  (LASA)  has  been  deployed  over 
a  large  •  -tion  of  southeast  Montana.  The  extent  of  this  array  is  about  an  order  of 
magnitude  greater  than  that  of  any  previous  instrument  used  in  the  study  of  seismic 
signals.  The  vertical  seismometers  of  the  experimental  LASA  are  spread  over  a  region 
which  can  be  contained  in  a  circle  having  a  radius  of  approximately  100  km.  The  array 
consists  of  21  subarrays  each  containing  25  seismometers  spread  over  a  circular  region 
with  a  radius  of  3.  5  km. 

Figure  1  shows  the  placement  of  the  21  subarrays  of  the  experimental  LASA. 

The  basic  configuration  of  seismometers  in  each  subarray  is  shown  in  Fig.  15B.  In 
many  subarrays  of  the  experimental  LASA  the  basic  subarray  geometry  was  rotated  up 
to  several  degrees.  Thus  subarray  configurations  are  essentially  identical  to  within 
rotations.  References  6  and  9  present  the  various  reasons  which  led  to  the  installation 
of  an  array  of  subarrays  and  to  the  particular  configuration  which  has  been  installed. 

Studies  are  currently  in  progress,  primarily  at  M.  I.T. ,  Lincoln  Laboratory, 
to  determine  how  best  to  utilize  the  voluminous  data  from  the  LASA  for  the  detection  of 
teleseismic  events  and  the  discrimination  between  natural  and  man-induced  tremors. 
Commensurate  with  these  goa.  s  the  array  is  also  being  used  to  study  aspects  of  earth 
structure,  propagation  phenomena,  and  seismic  noise  properties.  All  of  this  research 
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deals  with  the  use  of  data  obtained  from  an  already  existing  array.  The  prime  goal  is 
to  determine  the  ability  of  a  LASA  to  monitor  underground  nuclear  tests. 

If  it  is  determined  that  a  LASA  is  sufficiently  effective  for  seismic  detection 
and  discrimination  to  warrant  installation  of  several  more  such  arrays,  it  may  prove 
desirable  for  various  reasons  to  consider  geometries  different  from  that  of  the  experi 
mental  LASA.  It  is  with  this  contingency  in  mind  that  this  report  has  been  prepared. 
Several  families  of  possible  array  geometries  have  been  examined  to  find  a  few  of  the 
most  promising  ones.  Some  measures  of  the  sensitivity  of  array  properties  with  seis 
mometer  positions  have  also  been  obtained. 

The  same  computer  programs  used  to  study  and  evaluate  array  geometries  for 
this  report  could  also  be.used  to  evaluate  alternative  configurations  which  might  be 
suggested.  Also,  these  programs  might  well  be  useful  before  installation  to  study 
possible  detrimental  or  beneficial  effects  resulting  from  perturbations  of  an  intended 
array  geometry.  Such  perturbations  could  result  from  constraints  imposed  by  local 
topological  conditions,  land  availability,  or  other  factors. 
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II.  REVIEW  OF  BASIC  FACTS  ABOUT  ARRAY  PATTERNS 

The  effectiveness  of  an  array  and  associated  data  processing  scheme  for  the 
reduction  of  seismic  data  can  be  measured  by  the  response,  in  the  three  -dimensional 
frequency  wave  number  space,  of  the  array  and  associated  data  reduction  scheme. 
References  3,  4,  5,  6,  8,  and  10  are  six  of  many  in  which  the  representation  of  seismic 
signals  in  frequency  wave  number  space  is  discussed  and  the  interpretation  of  data  pro¬ 
cessing  methods  as  filters  in  that  space  is  explained. 

The  study  of  array  patterns  is  the  study  of  filter  characteristics  in  frequency 
wave  number  space  for  a  most  simple  kind  of  data  processing.  Specifically,  if  the 
array  contains  N  seismometers,  then  a  composite  trace  can  be  formed  simply  by  adding 
the  N  traces  and  dividing  by  N.  The  pattern  of  the  array  from  which  the  N  traces  are 
obtained  is  the  magnitude  of  the  gain  of  the  array  in  frequency  wave  number  space  when 
such  a  straight  sum  processing  technique  is  used. 

If  seismometer  responses  are  frequency  independent  and  no  frequency  dependent 
gains  are  inserted  before  summing,  then  the  gain  of  the  processing  scheme  is  independent 
of  frequency.  Under  these  conditions  only  the  response  in  two-dimensional  wave  number 
space  need  be  considered.  If  a  single  frequency  filter  operates  upon  the  sum  waveform 
or  identical  filters  operate  upon  every  one  of  the  individual  traces,  then  it  is  still  valu¬ 
able  to  consider  only  response  in  wave  number  space  assuming  no  filtering  before 
summing.  The  gain  at  any  point  in  frequency  wave  number  space  is  the  gain  in  the  plane 
where  frequency  is  zerc,  multiplied  by  a  function  of  frequency  only.  The  response  in 
the  plane  where  frequency  is  zero  is  the  response  in  two-dimensional  wave  number  space. 
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This  report  deals  only  with  the  gain  in  wave  number  space  obtained  when  straight  sum 
processing  is  employed.  Complex  processing  schemes  which  require  explicit  consid¬ 
eration  of  frequency  as  well  as  wave  number  in  order  to  be  characterized  are  beyond 
the  scope  of  this  paper. 

Different  patterns  can  be  obtained  as  a  function  of  array  geometry  without 
changing  the  data  processing  method.  A  comparison  of  the  array  geometries  can  then 
be  made  on  the  basis  of  the  patterns  which  represent  the  signal  processing  ability  of 
the  array.  Very  roughly,  one  wishes  to  synthesize  patterns  which  have  a  maximum 
gain  at  the  origin  in  wave  number  space  and  which  are  uniformly  low  in  some  region 
around  the  origin  of  wave  number  space.  The  extent  of  the  region  over  which  the  array 
must  have  low  gain  depends  upon  the  characteristics  of  seismic  noise. 

The  pattern  as  described  above  is  in  fact  the  magnitude  of  a  two-dimensional 
Fourier  transform  of  the  array  geometry.  Thus  the  pattern  of  an  array  with  vertical 
seismometers  at  positions  { x j , y j  :  j  =  1,...,N}  is 

1  N 

|H(kr,k  )|  =  jr  |  Yj  exp  [-12tt  (k  x  +k  y  )]  |  (1) 

a  y  j_|  a  j  y  j 

where  kx  and  k^  are  respectively  east  and  north  components  cf  horizontal  wave  number. 
H(~*  x,ky)  can  be  interpreted  as  the  two-dimensional  transform  of  an  impulsive  function 
with  unity  weights  on  each  impulse.  The  sifting  property  of  the  impulse  function  changes 
the  integrals  expressing  the  transform  into  the  above  sum  over  seismometer  positions. 

As  mentioned  in  the  introduction,  the  experimental  LASA  in  Montana  is  an  array 
of  roughly  identical  subarrays.  The  pattern  of  such  an  array  can  be  approximated  by 
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the  pattern  of  a  compound  array.  A  compound  array  is  one  which  is  an  array  of  identical 
subarrays.  A  compound  array  has  the  property  that  the  pattern  can  be  expressed  as  the 
product  of  the  pattern  of  a  subarray  and  the  pattern  of  a  configuration  having  a  single 


seismometer  located  in  the  place  of  the  subarrays.  To  se<_  this,  let  x  and  y  be  the 

pq  pq 

x  anu  y  coordinates  of  the  p^1  seismometer  in  the  q**1  subarray.  L,et 


X  =  X  +  X 

pq  p  q 


and 


o  , 

y  =  y  +y 
pq  p  q 


where  x°,y°  locate  the  center  of  the  pth  subarray  and  x  ,y^  locate  the  seismometers  cf 

p  p  q  q 

the  subarrays  with  respect  to  their  centers.  Let  P  be  the  number  of  subarrays  and  Q 
be  the  number  of  elements  in  the  subarray.  Then 


H(k  ,k  ) 
x  y 


_1_ 

PQ 


p  Q 

Z  Z  exp{-i2n[k  (x°+x  )  +  k  (y°+y  )]} 
p=J  q=l  x  p  q  y  P  q 


exp  [  -i  2tt  (k 


o 

x 

X  p 


+  k  y°)  ] 

y  p 


exp  [  -i  2tt  (k  x 
r  xq 


+  k  y  )] 

y  q 


Clearly  the  magnitude  of  the  product  of  the  right  is  the  product  of  the  magnitudes  and 
the  two  terms  of  the  product  are  the  patterns  of  the  array  centers  and  the  subarray 
respectively.  Henceforth  we  shall  use  the  terms  full  array  pattern,  array  pattern,  and 
suharray  pattern.  Full  array  pattern  and  array  pattern  are  used  to  distinguish  between 
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the  pattex  ->  generated  using  all  seismometers  in  an  array  and  the  pattern  generated  by 
the  centers  of  the  subarrays  in  a  compound  array,  respectively.  Subarray  pattern  of 
course  refers  to  the  pattern  of  a  subarray  in  a  compound  array. 

In  this  report  a  direct  confrontation  of  the  full  array  problem  has  been  avoided. 
Only  idealized  compound  arrays  having  the  factorization  property  just  described  have 
been  considered.  Subject  to  any  minor  positioning  errors  which  occur  during  installa¬ 
tion  and  even  to  subarray  rotations  such  as  those  or  the  experimental  LASA  this  separate 
consideration  of  array  and  subarray  patterns  rather  than  full  array  patterns  has  been 
found  to  be  a  valid  procedure,  as  will  now  be  demonstrated. 

Figure  2  shows  the  array  pattern,  evaluated  along  a  line  in  wave  number  space, 
of  the  actual  subarray  positions  of  the  experimental  LASA.  Figure  3  shows  a  similar 
plot  of  the  full  array  pattern  which  has  been  obtained  using  all  525  seismometers.  Note 
that  the  patterns  are  very  similar  for  wave  numbers  less  than  .05  cycles  Am  and  that 
for  larger  wave  numbers  the  full  array  pattern  has  smaller  sidelobes  than  the  array 
pattern.  The  location  of  subarrays  in  the  experimental  iLASA  was  fixed  assuming  that 
the  subarray  pattern  would  multiply  the  array  pattern  to  yield  the  full  array  pattern. 

The  idea  was  to  pick  a  subarray  and  array  so  that  this  product  would  yield  a  full  array 
pattern  having  all  but  the  nearest  lobes  of  the  array  pattern  significantly  reduced. 
Figures  2  and  3  clearly  show  that  such  a  philosophy  can  yield  a  satisfactory  full  array 
pattern  even  when  individual  subarrays  differ  from  each  other  by  being  rotated  several 
degrees. 
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The  experimental  LASA  cannot  have  a  pattern  which  is  exactly  the  product  of 
the  array  and  subarray  patterns  since  the  different  subarrays  have  patterns  differing 
by  rotations.  However,  such  an  approximation,  using  the  pattern  of  one  subarray,  is 
quite  good  for  wave  numbers  in  the  main  lobe  of  the  subarray  pattern  since  the  main 
lobe  changes  very  little  under  rotation.  For  larger  wave  numbers  the  approximation 
breaks  down  but  the  full  array  pattern  still  has  the  required  low  sidelobes.  F  gure  17B, 
if  the  wave  number  scale  is  multiplied  by  200/7,  is  a  cross  section  of  one  of  the  sub¬ 
array  patterns  for  the  experimental  LASA.  The  plots  in  Figures  2  and  3  extend  only 
about  as  far  in  wave  number  as  the  main  beam  of  a  subarray. 

Straight  sum  patterns  described  above  will  have  a  main  lobe  located  at  the  origin 
in  wave  number  space.  Thus  the  filter  operation  will  pass  signals  which  have  vertical 
incidence  upon  the  array  and  thus  have  infinite  horizontal  phase  velocity,  or  equivalent  y , 
zero  wave  number.  Other  signals  having  other  velocities  and  wave  numbers  will  be 
attenuated  as  indicated  by  the  pattern.  By  the  introduction  of  delays  into  seismometer 

outputs  before  summing  the  array  can  be  steered  to  pass  unattenuated  plane  waves 
arriving  from  any  specified  azimuth  having  any  specified  horizontal  phase  velocity.  At 
any  given  frequency  this  processing  has  a  gain  characteristic  in  wave  number  space  with 
exactly  the  same  shape  as  the  straight  sum  pattern  but  which  is  shifted  to  have  its  peak  at 
the  location  of  the  signal  at  that  frequency.  Thus  the  straight  sum  pat'  ;rn  is  a  useful 
measure  of  the  selectivity  of  any  array  even  when  the  array  is  steered  to  a  signal  having 
other  than  infinite  phase  velocity. 
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m.  PATTERN  REQUIREMENTS  AND  APPROACHES  TO  SYNTHESIS 


The  acceptability  of  patterns  depends  greatly  upon  the  noise  which  is  to  be 
rejected  when  me  array  is  sheered  to  an  event  of  interest  and  upon  the  expected  location 
in  frequency  wave  number  space  of  the  interesting  events.  Figure  4  is  a  map  of  the 
world  as  seen  by  an  array  centered  on  Miles  City,  Montana,  when  only  the  initial  P- 
phase  signal  at  one  Hertz  is  considered.^  Since  the  initial  P-phase  arrival  usually 
lias  most  frequency  content  at  about  one  Hertz  that  map  is  a  good  picture  of  the  location 
of  teleseismic  signals  of  interest  and  of  noise  which  results  from  distant  seismic  ac¬ 
tivity  and  reaches  the  array  via  earth  paths  similar  to  those  of  interesting  signals. 
From  the  consideration  of  such  a  picture  of  the  ^arth  as  seen  in  a  slice  through  fre¬ 
quency  wave  number  space  it  has  been  decided  that  any  LASA  which  is  to  be  used  for 
the  detection  of  underground  nuc  ear  tests  should  have  a  main  beam  with  a  radius  of 
about  .005  cycles  Am  or  less  in  wave  number  space.  That  is,  if  the  gain  of  the  LASA 
is  measured  at  a  distance  of  more  than  .005  cyclesAm  from  the  point  in  wave  number 
space  at  winch  it  is  steered,  then  the  gain  must  be  significantly  lower  than  unity.  This 
appears  to  allow  for  a  sufficient  degree  of  selectivity  in  steering  the  beam  to  some 
point  while  not  passing  similar  signals  from  other  parts  of  the  globe. 

If  subarray  or  array  geomttries  become  very  regular,  one  must  consider  a 
phenomenon  known  as  aliasing.  If  a  high  frequency  time  function  is  sampled  at  too  low 

a  rate,  the  sampled  waveform  can  appear  exactly  the  same  as  the  sampled  waveform 

2 

for  a  much  lower  frequency  signal.  Similar  problems  arise  when  spatial  signals  are 

3 

sampled  on  a  regular  grid.  Then  the  confusion  is  in  wave  number.  For  example, 


Figure  17B  is  a  section  of  the  pattern  of  the  regular  geometry  shown  in  Fig.  15B. 

Periodic  signals  coming  from  the  east  with  wave  numbers  greater  than  .  04  cycles  Am 
cannot  be  distinguished,  after  spatial  sampling  by  the  array,  from  some  other  signals 
having  wave  number in  the  range  0  -  .  04  cycles  Am.  TTiere  is  a  bad  sidelobe  at 
.08  cycles  Am,  exactly  twice  the  wave  number  at  which  aliasing  first  becomes  a  prob¬ 
lem.  Similar  lobes  appear  at  all  multiples  of  .08  cycles  Am  in  wave  number  space. 

Even  small  changes  in  the  spatial  sampling  positions  can  eliminate  true  aliasing 
as  a  real  problem.  However,  severe  sidelobe s,  almost  as  large  rs  unity,  can  still 
exist  in  regions  too  close  to  the  main  beam  of  the  pattern.  In  this  report  attention  has 
been  focussed  upon  such  lobes,  known  as  grating  lobes,  rather  than  upon  aliasing.  This 
was  done  since  it  is  the  grating  lobes  which  can  cause  trouble,  not  aliasing,  when  arrays 
are  not  extremely  regular. 

Consideration  of  possible  and  expected  noise  sources  has  led  to  the  conclusion 

that  a  full  array  pattern  should  have  no  grating,  lobes  closer  than  one  to  two  cycles  Am 
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to  the  center  of  the  main  beam  in  wave  number  space.  *  The  problems  of  aliasing 
are  avoided  by  assuming  a  sufficiently  irregular  placement  of  subarrays  and  of  seis¬ 
mometers  within  a  subarray. 

A  variety  of  approaches  to  the  problem  of  synthesizing  satisfactory  patterns 
using  about  500  seismometers  do  exist.  They  all  have  in  common  the  characteristic 
that  they  all  use  some  kind  of  trial  and  error  method  combined  with  reasonable  intuition. 
One  approach,  the  one  which  we  have  used  most  extensively,  is  to  restrict  the  arrays 
to  be  compound  arrays  similar  to  the  experimental  LASA.  The  full  array  pattern  can 
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be  synthesized  by  almost  independent  synthesis  of  the  array  and  subarray  patterns.  In 
particular,  the  array  pattern  might  be  synthesized  using  21  points  to  yield  a  .005  cycles 
per  km  half  beamwidth  or  less,  and  to  have  no  bad  sidelobes  closer  than  .  1  cycles /km 
to  the  center  of  the  main  beam.  This  choice  of  .  1  cycles /km  appears  to  give  a  reason¬ 
able  tradeoff  between  requirements  upon  the  array  and  the  subarrays.  The  subarray 
can  then  be  designed  to  have  a  main  beam  radiu~  of  .  1  cycles /km  or  less  and  no  bad 
sidelobes  closer  than  one  to  two  cycles/km.  Such  a  compound  array  would  then  meet 
the  design  requirements.  The  design  of  array  or  subarray  patterns  using  from  20  to 
25  points  must  proceed  via  trial  and  evaluation  methods.  Relatively  algorithmic  im¬ 
provement  methods  such  as  are  presented  in  Section  VII  can  at  best  be  of  value  to  find 
perturbations  of  a  geometry  which  would  yield  an  improved  pattern.  One  is  stymied 
when  no  nearby  geometry  will  yield  a  better  pattern  but  there  may  be  a  very  different 
geometry  having  a  much  more  satisfactory  pattern.  Tradeoffs  between  the  require¬ 
ments  imposed  upon  the  array  and  subarray  patterns  are  discussed  in  Section  VI. 

One  might  attempt  full  array  pattern  synthesis  using  some  of  the  formalized 
trial  and  error  methods  which  have  been  developed  for  the  synthesis  of  radar  and  sonar 
antenna  patterns.*  The  synthesis  of  a  LASA  geometry  is  very  similar  to  the  synthesis 
of  an  extremely  thinned  array.  However,  given  the  number  of  seismometers,  the  width 
of  the  main  beam,  and  constraints  upon  the  location  of  the  l.rst  grating  lobes,  the 
requirements  of  the  LASA  are  very  much  different  from  those  of  radar  arrays.  In  the 
radar  case  sidelobe  levels  less  than  20  db  down  from  the  main  beam  would  be  considered 
excessive.  Such  low  sidelobe  levels  can  be  maintained  for  thinned  arrays  only  by 
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increasing  the  width  of  the  main  beam.  For  the  LASA  such  an  increase  in  main  beam 
width  is  not  allowed.  Neither  is  a  significant  increase  in  the  number  of  seismometers. 
The  L.ASA  can  then  satisfy  the  design  criteria  only  by  a  significant  increase  in  sidelobe 
level.  It  is  thus  our  opinion  that  application  of  the  methods  suggested  for  density 
tapering  of  radar  arrays  cannot  yield  sidelobe  levels  significantly  lower  than  those 
which  have  been  achieved  using  compound  arrays.  Since  direct  synthesis  using  all  of 
the  seismometers  would  not  constrain  the  array  to  be  a  compound  array,,  such  an 
approach  has  not  been  pursued.  The  advantages  of  compound  arrays  (see  Section  IV) 
are  sufficiently  great  to  reject  a  direct  synthesis  of  a  full  array  since  the  direct  synthe¬ 
sis  doe 2  not  promise  any  significant  reduction  of  sidelobe  level  in  the  region  of  wave 
number  space  of  interest. 
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IV. 


COMPOUND  ARRAYS:  FAMILIES  OF  ARRAY  GEOMETRIES 


Many  engineering  advantages  are  gained  by  the  use  of  compound  arrays. 
Considerations  related  to  the  design  of  LASA  electronics  modules,  communication 
channels,  and  power  distribution  systems  all  argue  strongly  in  favor  of  the  use  of  com¬ 
pound  arrays.  The  use  of  such  arrays  also  provides  natural  divisions  for  data  pro¬ 
cessing  schemes.  Individual  subarrays  can  be  processed  and  then  the  results  combined 
in  a  second  stage  of  processing.  Failures  in  a  single  subarray  need  not  effect  the  re¬ 
duction  of  data  from  another  subarray. 

For  convenience  in  establishing  power  systems  and  communication  channels  it 
is  desirable  that  both  array  and  subarrays  be  set  out  along  several  radial  arms  connected 
to  a  central  point.  Such  a  configuration  tends  to  keep  the  length  of  lines  or  number  of 
separate  communication  channels  within  bounds.  The  use  of  configurations  which 
complicate  the  power  distribution  or  communication  networks  should  be  used  only  if 
considerable  increase  in  the  capabilities  of  the  array  can  be  gained.  No  such  geome¬ 
tries  have  yet  been  discovered. 

All  other  considerations  being  equal ,  arrays  graded  to  be  more  dense  towards 
the  center  are  desired  over  others.  This  grading  is  a  hedge  against  possible  loss  of 
coherence  of  teleseismic  P-phase  arrivals  as  the  signal  propagates  across  the  array. 

Array  geometries  having  the  most  remote  subarrays  on  a  circle  of  radius 
100  km  have  been  ctudied.  Such  a  dimension  was  picked  since  an  array  of  approximately 
that  size  must  be  used  in  order  to  obtain  the  narrow  main  beam  which  is  required.  An 
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array  should  have  at  least  three  subarrays  out  at  the  100  km  radius  circle.  More  than 
the  minumum  three  are  desirable  to  allow  the  use  of  redundancy  in  the  determination  of 
epicenters  using  a  limited  network  composed  of  a  single  output  from  each  of  the  sub- 
arrays  on  the  100  km  circle. 

Subject  to  the  above  constraints  a  relatively  extensive  search  for  possible 
suitable  locations  of  subarrays  has  been  completed.  Figures  5  through  8  show  typical 
members  of  the  families  investigated  as  possible  deployments  for  21  subarrays.  The 
number  of  points  on  a  ring  and  the  number  of  rings  is  fixed  for  each  family.  The 
spacing  of  rings  is  geometric  with  the  ratio  of  the  radius  of  adjacent  rings  being  a  con¬ 
stant.  This  constant  has  been  varied  to  generate  the  members  of  each  family.  The 
same  configurations,  but  with  arithmetic  spacing  of  rings,  have  also  been  evaluated. 

A  few  arrangements  with  subarrays  on  spiral  rather  than  straight  arms  and  with  ran¬ 
dom  placement  of  subarrays  were  considered  and  discarded  on  the  basis  of  their 
patterns.  Except  for  a  few  trials  the  number  of  subarrays  has  been  maintained  at  21, 
the  number  in  the  experimental  LASA  which  has  been  installed  in  Montana.  This  appears 
to  allow  a  reasonable  tradeoff  between  the  number  of  subarrays  and  the  number  of  seis¬ 
mometers  per  subarray.  There  is  no  reason  at  the  present  time  to  believe  that  another 
number  of  subarrays,  keeping  the  total  number  of  seismometers  constant,  could  yield 
significantly  better  performance  by  a  LASA. 

la~ge  amount  of  effort  has  l)een  expended  to  obtain  subarray  geometries  as 
opposed  to  array  geometries.  The  question  of  subarray  geometry  is  discussed  in 
Section  VI.  Briefly,  one  of  the  array  geometries,  suitably  scaled,  could  be  used  for  a 


subarray  with  21  seismometers,  or  other  subarray  geometries  using  a  different  number 
of  seismometers  could  be  employed. 
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V. 


BEST  ARRAY  GEOMETRIES  AND  THEIR  PATTERNS 


Figures  9  through  14  show  the  geometries  and  associated  patterns  of  the  most 
promising  arrays  found  in  each  of  the  families  of  arrays  described  in  the  preceding 
section.  The  families  with  five  or  seven  subarrays  on  each  ring  appear  to  be  the  most 
desi1  le. 

Array  patterns  have  been  graphically  presented  in  two  forms.  The  pattern  is  a 
gain  function  defined  on  the  two-dimensional  space  of  horizontal  wave  number.  Such 
patterns  have  been  graphically  presented  in  figures  such  as  9B  as  contour  plots.  The 
lines  on  the  plot  represent  the  level  of  the  pattern.  These  plots  were  obtained  by 
smoothly  connecting  the  points  on  a  plot  printed  by  the  line  printer  attached  to  a  large 
scale  digital  computer.  A  listing  of  a  Fortran  II  program  which  can  be  used  to  generate 
the  printer  plots  and  a  description  of  its  use  and  operation  are  contained  in  Appendix  A. 

Plots  such  as  that  shown  in  Figure  9C  were  obtained  by  evaluating  the  pattern 
along  a  ray  leaving  the  origin  of  wave  number  space.  The  plots  are  thus  cross  sections 
through  the  array  patterns.  Restriction  to  only  a  few  radial  directions  and  the  use  of 
a  more  accurate  output  device  than  a  printer  plot  allowed  a  mere  accurate  and  detailed 
picture  of  the  pattern  to  be  obtained  in  directions  of  particular  interest.  Plots  were  ob¬ 
tained  in  tho°e  directions  in  which  the  cruder  printer  plots  indicated  the  most  severe 
sidelobes  v  .d  be  found.  A  listing  of  Fortran  IV  programs  used  to  obtain  cross  sec¬ 
tional  plots  is  given  in  Appendix  B.  That  appendix  also  describes  the  operation  and  use 
of  those  programs.  The  programs  have  been  constructed  so  that  minimal  modifications 
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would  allow  one  to  obtain  cross  sectional  patterns  without  the  use  of  a  Stromberg 
Carlson  4020  plotter  which  was  the  output  device  used. 

The  programs  described  above  can  also  be  easily  modified  to  help  interpret  in 
wave  number  or  frequency  wave  number  space  the  effect  of  processing  schemes  much 
more  complex  than  the  straight  sum  processing  used  for  the  evaluation  of  alternative 
geometries.  As  mentioned  in  Section  II  the  interpretation  of  more  complex  processing 
in  frequency  wave  number  space  is  outside  the  scope  of  this  report  and  has  not  been 
included. 

In  each  family  the  arrays  were  judged  most  satisfactory  when  they  had  the  least 
severe  sidelobes  in  the  range  from  zero  to  .  1  cycles /km  and  had  beam  widths  which 
were  sufficiently  small.  As  discussed  above  the  beam  width  should  be  .  01  cycles  Am 
or  less.  Sidelobes  further  than  .  1  cycles  /km  from  the  main  beam  will  be  sufficiently 
attenuated  by  subarrays  if  the  subarrays  are  properly  set  out.  Since  no  single  scalar 
quantity  can  adequately  represent  the  quality  of  the  various  arrays,  the  choice  of  best 
was  necessarily  somewhat  subjective.  At  *  .at  one  representative  best  member  from 
each  major  family  imestigated  has  been  saved. 

Comparisons  between  members  of  the  same  family  or  of  different  families  of 
array  geometries  is  made  a  bit  less  complex  by  noting  certain  symmetries  in  patterns 
and  invariances  with  geometries.  Ncte  that  contour  maps  of  patterns  have  been  drawn 
only  over  relatively  restricted  pie  shaped  sections  of  wave  number  space.  Only  those 
sections  are  necessary  since  the  entire  pattern  is  obtainable  simply  by  completing  the 


circle  with  identical  pie  slices.  This  result  and  certain  invariances  for  patterns  of 
arrays  with  odd  numbers  of  subirrays  per  ring  are  described  in  more  detail  in 
Appendix  C. 

The  array  configuration  shown  in  Figure  13A,  which  has  been  picked  as  the 
best  geometry  with  four  subarrays  per  ring,  yields  a  pattern  which  appears  to  be  defin¬ 
itely  inferior  to  those  which  use  three,  five  or  seven  subarrays  per  ring.  The  pattern 
for  the  geometry  shown  in  Figure  13A  has  significant  sidelobes  which  are  no  further 
than  .015  cycles  Am  away  from  the  center  of  the  main  beam.  Even  more  severe 
sidelobes  appear  at  .05  and  .07  cycles/km  from  the  main  beam.  This  best  member  of 
the  family  of  arrays  with  four  subarrays  per  ring  is  quite  similar  to  the  actual  config¬ 
uration  of  the  experimental  LASA  which  has  been  installed  in  Montana.  Similar 
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severe  sidelobes  have  been  noted  for  that  experimental  array.  *  *  These  can  be 
seen  in  Figure  2. 

It  is  our  opinion  that  the  array  shown  in  Figure  14A  which  has  three  subarrays 
per  ring  can  also  be  rejected  as  the  most  valuable  array.  The  pattern  for  this  array 
has  a  broader  main  beam  than  that  of  other  geometries  under  consideration.  Sidelobes 
are  not  excessive,  but  are  slightly  larger  in  amplitude  than  those  of  alternative  patterns. 
If  the  array  is  stretched  over  a  larger  area  to  reduce  the  width  of  the  main  beam,  then 
the  sidelobes  will  move  in  closer.  Finally,  more  than  three  subarrays  are  desirable 
on  the  largest  ring  for  the  purpose  of  determining  epicenters  using  the  array. 

The  configuration  shown  in  Figure  10A  is  probably  to  be  preferred  to  that  of 
Figure  9A  which  also  has  seven  subarrays  per  ring.  The  array  in  Figure  9A  has 
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subarray,?  concentrated  more  towards  the  100  km  radius  ring.  Since  the  validity  of 
patterns  depends  somewhat  upon  the  coherence  of  a  signal  across  the  array,  *+iis  more 
spread  out  array  would  probably  be  more  adversely  affected  by  loss  of  coherence. 
Indeed ,  any  signal  processing  scheme  which  depends  upon  signal  coherence  will  be 
more  s3iiously  affected  by  loss  of  coherence  when  subarrays  are  spread  out  than  when 
they  are  clustered  together.  One  might  consider  scaling  down  the  size  of  the  array 
shown  in  Fig.  9A  since  it  appears  tc  have  a  very  narrow  main  beam.  This  should  be 
avoided  since  it  will  seriously  effect  the  usefulness  of  the  outer  ring  of  the  array  for 
the  determination  of  epicenters.  The  decision  between  these  two  arrays  is  at  best 
subjective. 

A  compari  on  of  the  geometries  of  Figures  11A  and  12A,  having  five  subarrays 
per  ring,  with  the  configurations  having  seven  per  ring  is  completed  in  Section  VII. 
That  comparison  is  deferred  because  the  application  of  a  gradient  technique  described 
there  has  resulted  in  some  changes  in  the  sidelobes  of  the  pentagonal  arrays. 


18 


VI.  TRADEOFF  BETWEEN  SUBARRAYS  AND  ARRAY  IN  A  COMPOUND  ARRAY 

In  the  above  sections  it  has  been  rather  casually  assumed  that  any  bad  sidelobes 
of  the  array  pattern  farther  out  than  .  1  cyclesAm  could  be  sufficiently  attenuated  by 
the  subarrays  to  be  negligible.  It  has  even  been  implicitly  assumed  that  some  attenu¬ 
ation,  say  a  few  db,  might  be  obtained  at  wave  numbers  having  magnitude  as  small  as 
.05  cyclesAm  or  less.  These  assumptions  are  discussed  in  this  section  and  the  im¬ 
portance  of  the  relative  sizes  of  array  and  subarrays  is  emphasized. 

Assume  tu?t  the  array  geometry  shown  in  Fig.  10A  is  used  and  that  each  subarray 
has  seismometers  in  the  same  configuration  but  that  the  largest  ring  of  seven  seismom  - 
eters  has  a  10  km  diameter.  The  subarray  is  then  1/20  the  size  of  the  array.  The  gain 
of  such  a  compound  array  is  found  by  forming  the  product  of  the  array  and  subarray 
gains. 

Consider  the  most  severe  array  sidelobe  shown  in  the  center  frame  of  Fig.  10C. 
That  lobe  is  about  .037  cycles/km  from  the  origin  and  has  level  .42.  The  subarray 
pattern  will  have  the  same  shape  but  the  wave  number  scale  will  be  multiplied  by  the 
ratio  of  array  to  subarray  diameters.  If  the  subarray  has  diameter  10  km,  it  will  be 
able  to  reduce  the  above  mentioned  array  lobe  to  a  level  of  about  .  3b.  The  use  of  a 
subarray  having  a  diameter  of  15  km  would  attenuate  the  lobe  even  more  to  abo”  .  3. 

If  the  subarray  diameter  were  reduced  to  7  km,  this  same  lobe  would  remain  larger 
than  about  .38.  These  figures  are  approximate  only  because  subarray  gains  have  been 
obtained  from  Fig.  10C  by  changing  the  wave  number  scale.  The  point  is  that  even  for 
wave  numbers  as  small  as  .  037  cycles/lm  some  meaningful  attenuation  can  be  obtained 
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from  subr  ys  having  diameters  from  10  to  15  km.  The  attenuation  from  a  7  km 
subarray  is  very  slight  so  close  to  the  origin  of  wave  number  space. 


The  larger  subarray  diameters  maintain  similar  advantages  all  the  way  out  to 
.  1  cycles  Am.  The  array  sidelobes  are  more  effectively  attenuated  by  subarrays  with 
diameters  10  to  15  km  than  by  those  with  diameter  7  km.  This  will  tend  to  be  true 
independently  of  the  particular  subarray  configuration  which  is  used. 

If  the  configuration  of  Fig.  10A,  scaled  to  have  a  diameter  between  10  and 
15  km,  is  used  as  a  subarray,  then  the  .42  level  closest  sidelobe  will  appear  between 
.48  and  .74  cycles  Am  from  the  center  of  the  subarray  pattern.  This  sidelobe  is  not 
extremely  detrimental  and  may  be  tolerable  in  the  subarray  pattern. 

Figure  17B  is  a  cross  section  of  the  pattern  generated  by  the  geometry  of  a 
subarray  of  the  experimental  LASA.  The  configuration  has  been  scaled  to  a  200  diam¬ 
eter  circle  for  convenient  comparison.  The  true  subarray  lies  w.  min  a  circle  with  a 
diameter  of  7  km.  Figure  15B  shows  the  seismometer  configuration.  Observe  that  the 
pattern  for  this  configuration  has  a  close  in  sidelobe  at  level  .  44,  higher  than  the  one 
discussed  above  for  a  subarray  using  only  21  elements.  Even  if  the  subarray  has 
diameter  7  km,  this  lobe  is  at  .4  cycles  Am.  In  general,  the  pattern  drawn  in  Fig. 

10C  is  definitely  superior  to  that  drawn  in  Fig.  17B.  The  superiority  of  the  21 -element 
subarray  in  terms  of  sidelobes  is  maintained  even  when  its  size  is  increased  to  10-15 
km  while  the  array  using  25  seismometers  is  maintained  at  7  km.  The  21 -element  sub- 
array  has  a  generally  lower  sidelobe  level  and  has  it  ™er  a  larger  range  of  wave  num  - 
bers  than  t».e  25  element  subarray  does. 
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’-T? 


It  seems  that  the  use  of  21  seismometers  in  a  subarray,  placed  on  three  rings, 
containing  seven  seismometers,  may  yield  a  subarray  with  characteristics  superior  to 
those  of  the  LASA  Montana  subarrays  which  use  25  seismometers.  One  might  then  ex¬ 
pect  that  a  subarray  configuration  which  uses  25  seismometers  and  is  definitely  superior 
to  either  of  those  discussed  above  could  be  found.  No  extensive  search  for  such  a 
geometry  has  been  completed.  However,  one  family  of  geometries,  having  five  lings 
spaced  geometrically  and  each  having  five  seismometers  per  ring,  has  been  briefly 
investigated.  Figure  15/i  shows  the  resulting  subarray  end  Figs.  16  and  17A  show 
pattern0  The  lobe  of  level  .49  at  a  distance  of  .05  cycles/km  from  the  center  of  the 
main  beam  is  the  worst  characteristic  of  the  pattern.  This  lobe  might  be  removed  using 
sensitivity  functions  as  described  in  the  next  section  or  a  completely  different  family 
using  25  seismometers  might  be  investigated.  The  synthesis  of  configuration  using 
25  points  has  not  been  pursued. 
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VH.  SENSITIVITY  AND  PATTERN  IMPROVEMENT 

The  sensitivity  function  of  a  pattern  or  array  is  a  measure  of  how  much  the 
pattern  will  be  effected  by  small  changes  in  seismometer  or  subarray  positions.  There 
are  two  principal  reasons  for  obtaining  a  sensitivity  function.  First,  of  course,  such 
a  functior  indicates  how  accurately  subarrays  must  be  located  so  that  the  pattern  will 
not  be  seriously  modified  by  errors  in  placement.  Such  a  functio  n  can  also  be  used  to 
help  find  improved  geometries.  It  can  indicate  how  the  geometry  should  be  perturbed 
to  achieve  certain  desired  pattern  changes.  A  useful  sensitivity  measure  is  defined  in 
the  following  paragraphs. 

Let  r^  be  the  two-dimensional  vector  which  locates  the  subarray  in  an  array 
or  qtn  seismometer  in  a  subarray.  Suppose  that  u^  is  a  unit  vector  which  indicates  the 

ft.  _ 

direction  in  which  the  q  element  might  be  moved  and  that  F(r^)  is  a  complex  valued 

function  of  r^.  Tne  function  F  can  of  course  be  a  pattern  with  dependence  upon  the 
th 

position  of  the  q  point  of  the  array  made  explicit.  If  it  exists,  the  following  directional 
derivative  can  be  defined. 


F(r  +eu  )  -  F(r  )  9F(r  +eu  ) 

#  =  um  —a — a - Lai  =  _la: — al 

*q  e  |  o  e  e  =  0 


It  is  a  matter  of  straightforward  substitution  to  verify  that  if 


F  =  F  F 
1  2 


22 


*  -F  5  +  f  5 

Bu  1  Bu  2  Bu 

q  q  q 


which  is  the  usual  rule  of  differentiation  of  products.  Other  rules  such  as 


BF  _,s  ~  1  BF 

aT  =  sF  sT 
q  q 


also  hold  for  the  directional  derivative. 


A  physical  meaning  can  easily  be  attached  to  BF/Bu^.  To  a  first  order 


approximation  we  have 


3F(r  ) 


F(r  +eu  )  =  F(r  )  +  e  —=■ 
q  q  q  Bu 


if  c  is  sufficiently  small.  If  F  is  a  pattern,  then  e  [BF(r^)/Su^]  is  roughly  the 
change  in  the  pattern  when  subarray  or  seismometer  q  is  moved  e  km  in  the  direction 


of  u  . 

q 


Let  (k,u^)  an<*  (k,r^)  denote  usual  scalar  products  of  the  indicated  vectors. 


Then  the  magnitude  of 


,  in  _ 

H  =  rr  E  exP  [  “i  2tt  <k,r  )  ]  , 
j=l  J 


which  equals  (HH*)1^,  is  the  pattern  of  an  array.  One  can  easily  obtain 
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From  this 


-i  ~  <k*uq)  exp  [-i2iT  <k,rq> 


-f=S“  =  “  (k»uq)  Im  {H  exp  [i  2tt  (k,rq>  ]  } 


W*-  =  “af2-  =  "  Fpf  <k’V  to{He*Pli2"<k-rq>l} 


Considered  as  a  function  of  k  this  is  a  sensitivity  function  which  reflects  the  effects  of 

th  — 

moving  the  qr  point  on  an  array  in  the  u  direction. 

b!  H I 

It  is  convenient  and  useful  to  express  — W-L  in  terms  of  polar  coordinates  of 

duq 

all  quantities.  Let  |h|  and  cp  be  the  magnitude  and  phase  of  the  complex  gain  H.  Let 
0  be  the  angle  of  u  and  let  r  ,  tlf  be  the  magnitude  and  direction  of  the  position  r  . 

q  q  q  q  q 

_  3 1 H  j 

Finally  let  k  and  0  be  magnitude  and  direction  of  k.  Using  this  notation,  L  -»■  can  be 

uq 

expressed  as 


^  k  cos  (0-0q)sin  [cp  +  2nkrqcos  (0-t|fq)]  • 


It  is  this  expression  using  polar  notation  which  has  been  most  useful. 

It  is  now  very  easy  to  obtain  an  upper  bound  npon  the  possible  effects  of  moving 


a  seismometer  or  subarray.  In  particular , 
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2rrk 

N 


lU'L 

Su 

q 

This  bound  on  sensitivity  is  inversely  proportional  to  the  number  of  elements  and 

directly  proportional  to  the  size  of  the  wave  number  at  ‘‘he  point  in  wave  number  space 

which  is  of  interest.  Small  changes  in  subarray  positions  can  greatly  effect  the  pattern 

for  large  values  of  wave  numbers.  Note  also  that  this  bound  is  independent  of  r^,  the 

th 

magnitude  of  the  position  vector  for  the  q  n  element.  Roughly  speaking,  a  change  of 
e  km  in  the  position  of  an  array  point  cannot  change  the  pattern  in  wave  number  space 
more  than  2nke/N.  Thus  a  change  of  1  km  in  the  position  of  one  element  of  a  21-element 
array  of  subarrays  cannot  change  the  pattern  by  more  than  (.  3)  k.  For  k  less  than 
.  i  cycles  Am,  this  is  less  than  .  030.  Based  upon  this  bound  one  can  conclude  that 
changes  as  large  as  1  km  in  the  location  of  subarrays  may  not  result  in  drastic 
changes  in  patterns.  Certainly  the  pattern  of  a  modified  array  must  be  checked  but  it 
is  hoped  that  changes  as  large  as  one  kilometer  in  position  can  be  tolerated  without  un¬ 
due  concern.  Similar  computations  for  subarrays  covering  a  region  having  a  10  km 
diameter  indicate  that  errors  on  the  order  of  .  05  km  may  not  have  too  adverse  effects 
although  the  effects  must  be  checked.  Of  course,  if  these  bounds  are  divided  by  the 
number  of  elements,  then  the  check  is  not  needed  as  long  as  changes  as  large  as  .  03 
in  pattern  are  considered  acceptable.  Such  an  approach  requires  much  more  accurate 
placement  of  elements. 

It  may  not  be  physically  possible  to  locate  a  subarray  less  than  a  kilometer  away 
from  a  point  which  is  specified  for  an  idealized  array.  Such  a  situation  need  cause  no 
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serious  difficulty.  Suppose  some  realizable  alternative  to  an  idealized  array 
configuration  is  under  consideration.  The  pattern  for  the  configuration  can  be  quite 
completely  obtained  in  the  regions  of  wave  number  space  of  interest  using  programs 
such  as  are  listed  in  the  appendices.  A  study  including  two  contour  plots  and  a  dozen 
cross  sectional  plots  would  require  only  about  five  minutes  of  time  on  an  IBM  7094 
computer.  The  cost  of  such  a  study  of  an  array  configuration  is  so  slight  that  it  can 
always  be  done  for  any  configuration  under  serious  consideration. 

Since  it  is  so  inexpensive  to  evaluate  a  pattern,  it  seems  reasonable  to  place 
bounds  on  allowed  variations  in  subarray  locations  by  specifying  a  surface,  defined  in 
wave  number  space,  below  which  the  array  pattern  must  lie.  Such  a  requirement  could 
be  used  in  conjunction  with  quite  weak  requirements  upon  allowed  changes  in  subarray 
position  from  idealized  positions.  For  example,  changes  of  position  of  the  order  of 
5-10  km  would  not  be  unacceptable  if  the  pattern  was  not  deteriorated  by  those  changes. 
If  such  a  method  of  specification  is  used,  the  strong  requirements  on  position  imposed 
by  the  use  of  sensitivity  functions  can  then  be  relaxed. 

Some  success  has  been  achieved  in  using  sensitivity  functions  to  improve 
patterns  by  moving  subarray  locations.  A  Fortran  IV  program  has  been  written  which 
evaluates  sensitivity  of  patterns  with  respect  to  the  radius  of  a  ring  of  elements  in  an 
array.  This  is  done  simply  by  setting  9^  equal  to  \jr  and  summing  over  those  q  which 
index  the  elements  of  the  ring  of  interest.  The  program  computes  and  plots  this  sensi¬ 
tivity  function  for  points  along  a  ray  in  wave  number  space.  The  program  is  listed  and 
its  use  described  in  Appendix  D. 
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Figure  18  shov/s  the  sensitivity  of  a  pattern  to  changes  in  the  radii  of  the  five 
element  rings  of  the  array  shown  in  Fig.  12A.  The  sensitivity  has  been  evaluated  along 
a  direction  in  wave  number  space  in  which  are  found  the  most  severe  sidelobes  of  the 
pattern.  These  sensitivity  functions,  as  well  as  all  others  plotted  in  this  report, 
should  be  divided  by  the  number  of  elements  (21)  to  agree  with  sensitivity  as  defined 
above. 

It  was  determined  using  the  sensitivity  function,  that  the  bad  lobe  at  about 
.0425  cycles  Am  in  wave  number  could  be  reduced  by  reducing  the  radius  of  the  ring  of 
elements  having  a  25  km  radius,  judging  from  the  sensitivity,  this  change  should  re¬ 
duce  the  lobe  at  .07  cycles  Am  a  comparable  amount  if  not  more. 

Figure  19  shows  the  effect  upon  pattern  and  sensitivity  of  reducing  the  inner 
ring  radius  to  20  km.  The  plots  are  along  the  same  direction  in  wave  number  space  as 
those  of  Figure  18.  The  lobe  at  .0425  cycles  Am  has  been  reduced  from  .  52  to  .43, 
a  nontrivial  improvement.  The  second  severe  lobe  did  not  fare  sc  ;ell.  As  the  radius 
of  the  ring  decreased,  the  sensitivity  changed  in  such  a  way  that  the  net  effect  upon  the 
lobe  was  slight. 

Figure  20  shows  the  effects  upon  the  pattern  in  a  different  direction  in  wave 
number  space  of  the  changes  described  above.  The  direction  was  that  having  the  most 
severe  adverse  effects.  Specifically,  the  lobe  at  about  .055  cycles  Am  has  grown  to  a 
value  of  .  56.  This  degradation  was  roughly  predicted  by  the  sensitivity  function  eval¬ 
uated  along  the  direction  18°  north  of  east. 
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The  relative  quality  of  the  patterns  of  the  modified  and  unmodified  arrays  is  not 
clearly  defined.  It  is  our  opinion  that  the  reduction  of  the  level  of  the  nearest  sidelobe 
looking  eastward  is  worth  the  cost  paid  in  other  directions.  In  general  the  adverse 
effects  of  sidelobes  becomes  greater  as  they  get  closer  to  the  main  lobe.  Very  close 
sidelobes  cannot  be  killed  by  subarray  processing. 

No  other  modifications  cf  this  basic  geometry  have  been  attempted.  Sensitivity 
functions  have  not  indicated  any  changes  expected  to  yield  significant  pattern  improve¬ 
ment. 

A  procedure  similar  to  that  described  above  has  also  been  followed  starting 
with  the  array  shown  in  Fig.  11A.  The  results  are  summarized  in  Figs.  21,  22,,  and 
23.  By  increasing  the  ring  having  radius  55  km  to  57  km  the  lobe  at  about  .063  cycles 
per  kilometer  in  the  eastward  direction  has  been  reduced  from  .  5  to  .  43.  The  lobe  at 
about  .054  has  increased  from  .  35  to  .44.  The  relative  quality  of  the  modified  and 
unmodified  arrays  is  so  tenuous  that  no  attempt  is  made  to  rank  them. 

No  further  attempts  at  modifying  these  patterns  have  been  made.  It  is  possible, 
for  example,  by  decreasing  the  radius  of  the  30. 25  km  ring,  that  marginal  improve¬ 
ment  in  sidelobes  could  be  obtained.  The  small  amounts  of  expected  improvement  did 
not  seem  significant  enough  to  warrant  any  further  search. 

Figures  24  and  25  are  pattern  and  sensitivity  plots  for  the  arrays  of  Figs.  9A 
and  10A,  respectively.  The  directions  in  wave  number  space  were  picked  to  exhibit 
most  severe  sidelobe  activity.  No  changes  of  radii  of  rings  appear  to  promise  im¬ 
proved  sidelobes. 
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After  considering  ail  the  various  patterns  described  in  this  report,  it  is  our 
opinion  that  the  geometry  of  Tig.  10A  using  seven  subarrays  per  ring  is  the  best  to  use 
for  a  LASA  Unfortunately ,  it  does  not  appear  that  the  lobe  of  level  .  42  at  a  distance 
of  .  037  from  the  main  beam  can  be  removed  by  slight  changing  of  the  radii  of  rings. 
Even  with  that  lobe  this  geometry  seems  most  desirable.  No  other  bad  lobes  are 
present.  This  is  not  true  of  geometries  having  five  subarrays  per  ring.  Those  arrays 
have  first  significant  lobes  a  bit  farther  out  but  have  a  generally  higher  sidelobe  level 
and  have  other  very  significant  sidelobes  closer  than  .  1  cycles/km  from  the  center  of 


the  main  beam. 


VIII.  A  comparison  with  the  experimental  lasa 


Figures  2  and  IOC  show  cross  sections  of  the  patterns  for  the  experimental 
LASA  geometry  and  for  a  proposed  geometry  which  would  deploy  21  subarrays  on 
three  rings.  For  wave  numbers  in  the  range  .  02  -  .  2,  the  sideiobe  levels  are  roughly 
the  same  for  the  two  geometries.  However,  for  wave  numbers  less  than  .  02  in  mag¬ 
nitude  the  proposed  configuration  shown  in  Fig.  10C  yields  a  pattern  very  much  better 
than  the  pattern  of  the  experimental  LASA. 

The  experimental  LASA,  depending  upon  one's  point  of  vie\. ,  either  has  a  main 
beam  nearly  .  04  cycles  /km  in  diameter  or  has  a  bad  sideiobe  only  .015  cycles/km 
from  the  center  of  the  main  beam.  The  proposed  configuration  has  no  such  close  side- 
lobfc  and  has  a  main  beam  with  a  diameter  of  about  .  01  cycles/km  in  wave  number 
space.  The  pattern  for  the  proposed  geometry  shown  in  Fig.  10A  is  sufficiently  better 
than  the  pattern  for  the  experimental  LASA  to  warrant  serious  consideration  if  any  other 
LASAs  are  constructed.  It  should  be  noted  that  pattern  improvement  has  been  obtained 
primarily  in  the  important  region  of  wave  number  space  close  to  and  including  the  main 
beam. 
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IX. 


SUMMARY 


A  study  of  possible  geometries  for  a  Large  Aperture  Seismic  Array  (LASA) 
has  been  completed.  For  a  variety  of  reasons  the  array  was  restricted  to  being  com¬ 
posed  of  21  subarrays  deployed  over  a  circular  region  having  a  diameter  of  200  km. 
From  each  of  several  families  of  array  geometries  the  members  yielding  best  patterns 
were  saved  and  discussed.  Two  families,  those  having  three  or  four  subarrays  on 
each  of  several  rings,  were  discarded  as  definitely  inferior  to  other  families.  The 
relative  quality  of  arrays  using  five  or  seven  subarrays  per  ring  was  not  as  clear.  The 
geometry  shown  in  Figure  10A appears  to  be  the  best  all  around  compromise.  That 
geometry  yields  a  significantly  better  main  beam  than  the  experimental  LASA  does. 
Other  possible  geometries  and  their  patterns  have  been  saved  for  comparison. 

The  diameter  of  subarrays  in  the  experimental  LASA  in  Montana  is  7  km. 

Pattern  studies  indicate  that  this  should  be  increased  to  from  10  to  15  km  for  other 
installations.  Even  using  from  20  to  25  seismometers  per  subarray,  undesirable 
aliasing  and  unwanted  grating  lobes  can  be  avoided  by  the  use  of  suitable  subarray 
geometries. 

Sensitivity  functions  for  array  patterns  have  been  developed.  These  indicate 
that  position  changes  of  the  order  of  .  5%  of  the  diameter  of  an  array  or  less  should 
have  only  slight  effect  upon  the  array  pattern.  Fortunately  it  should  not  be  necessary 
to  maintain  such  accurate  placement  of  subarrays  with  respect  to  some  idealized  array. 
The  cost  of  evaluating  the  parern  of  a  configuration  is  small  enough  to  allow  this  to  be 
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done  for  any  configuration  under  serious  consideration.  Any  changes  in  a  configuration 
can  be  directly  evaluated  in  terms  of  the  pattern.  The  use  of  sensitivity  functions  to 
improve  geometries  has  been  investigated  with  moderate  success. 


A  PPENDIX  A 


PROGRAM -FOURIER  TRANSFORM  WITH  PRINTER  PLOT 


The  Fortran  II  program  listed  below  computes  the  pattern  defined  by  (1)  for 
values  of  k^k^  ^n  a  rectangular  grid  (100  x  100)  in  the  wave  number  plane.  Execution 
of  the  program  generates  a  printer  plot  which  is  a  map  of  the  patte  n.  If  instruction 
111  and  the  one  preceding  it  are  removed  from  the  program,  then  the  program  can  be 
used  to  generate  gain  plots  in  wave  number  space  resulting  from  weighted  sum  pro¬ 
cessing  rather  than  straight  sum  processing.  The  inclusion  of  these  instructions 
avoids  the  need  to  punch  a  one  on  each  card  which  locates  a  subarray. 

The  Plot 

The  plot  is  a  contour  map  over  a  rectangular  region  of  k  ,k  space.  One 

x  y 

corner  of  the  rectangle  is  at  position  (XK,  YK).  The  diagonally  opposite  corner  is  at 
(XK  +  RX,  YK  +  RY).  These  points  are  at  the  upper  left  and  lower  right  corners  of  the 
output  plot. 

The  user  specifies  a  parameter  DELTA  where  0. 05  <  DELTA  <  1.  The  range 

of  the  magnitude  of  gains  is  divided  into  intervals  [  (h- 1)  DELTA,  hDELTA]  with 

h  =  1,2,. . .  Those  points  (k  ,k  )  on  the  grid  for  which 

x  y 


(h-1)  DELTA  <  |  Gain  j  <  hDELTA 


are  marked  with  the  hcn  character  from  the  list 


(1, 2,3,4, 5, 6, 7, 8.  9,  A ,B,C, D,  E,  F. G,H,I,J,K, L.M.N, blank)  . 


In  order  to  aid  the  user  to  read  the  contour  plot,  blanks  are  substituted  for  characters 
when  several  contiguous  characters  are  the  same.  Only  the  left  most  of  such  a  group 
of  characters  is  written  in  any  row. 

Use  of  the  Program  (Fortran  II) 

Following  the  *DATA  card,  there  must  be  four  cards  containing  octal  coding, 
format  (6012),  for  the  characters  in  the  list  given  above.  Following  these  four  cards 
there  can  be  any  number  of  sets  of  data.  Each  set  of  data  contains  a  card  with  format 
(215, 5 F 10. 5)  containing  N,  NCTRi,  DELTA,  RX,  RY,  XK,  YK.  The  variable  N  speci¬ 
fies  the  number  of  seismometers.  NCTR1  is  a  control  number.  If  NCTR1  =  3,  we  ob¬ 
tain  a  plot  of  the  magnitude  of  gain.  If  NCTR1  =  2,  then  a  listing  of  N  times  the  gain 
is  given.  If  NCTR1  =  1,  both  the  plot  and  the  listing  are  given,  The  variables 
DELTA,  RX,  RY,  XK,  YK  have  been  described  above.  Within  each  set  of  data,  this 

card  is  followed  by  N  cards  with  format  (42X,  3F10. 3).  They  contain  (y..x.,w.), 

j  j  J 

j  =  1, . . .  ,N  where  Yj»Xj  are  the  north  and  east  components  of  the  location  of  the  jth 
seismometer  or  subarray  and  w^  is  the  weight  attached  to  it  for  weighted  sum  pro¬ 
cessing.  For  pattern  studies  the  w  need  not  be  punched  and  are  ret  to  unity  by 
instruction  111. 

The  four  octal  coded  cards  appear  only  once,  immediately  after  *DATA,  regard¬ 
less  of  the  number  of  sets  of  data  processed. 

The  computation  of  a  plot  using  21  seismometers  requires  about  one  minute  or 
less  of  time  o”  an  IBM  7094  computer. 
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T*P  FOURIER  TRANSFORM  WITH  PRINTER  PLOT - LACOSS  9/22/6) 

MOp|Fm?_NOV.«*1965-™  DOES.  NOT  HtUUIK;  RtAQINO  OF  MtlQHl* 


ASSUMES  all  rfK«TS  TO  *t  uhllf 

BiHSir  OA  A>fcOGUY>6O0l*S‘xT&Li24Wil»C:&-ltLl24Qll*VALUI4£&> 

DIMENSION  RNUMt6opj.*KlNm00ltALPMAt24|  _ _ 

DIMENSION  MibUOi 

REAy  INPUT  TAPE  2.2^0<,.l  ALPHAI  J».J*1  .24) 


t 

C 


FORM  5  IN  AND- COS  TABLES 
SO  i0.fi*U©6  _ 


FL I  ■  I 

THETA*  13.14159266/  I  2‘JO.)«TFl  1-1,  I 
_ 10  SNT t>L(  1 ) *SlNF  UilLI.AJ _ _ 

DO  12  1*1201*2401 
SNTHU  I  )«-SNTBU  1-1260) 

oo  i*  1'itim 

14  CSToj.  <  l  )  *SNT«L  l  i  *buo  > 

DO  16  l  ■  1002  »2**Ul 
16  CSTOLI  U*S?»JLttLll-ieyu| 


W  NEW?  HD  A 

READ  Uif*y?_U££  i.lOVt»NjNCIBliOtkI4*J4Aft6E*1gA«6fclt*fcST  wTF.Sl 

1000  F0RHATf2I5.5F10.i > 


•RITE  OUTPUT  TAPt  J.Juww.N* 

DELTA  .RANOE  X  .RANGE  Y.XX6T*Y»V? 

3UU0  FORMAT  4?HlVAL5JtS  OF  N #ULL T  A .RANGtX  »K  ANGE.  Y  *  A&ST  » YLST/  ilO.5FlS.2i 

- 

c 

c 

N--NUKGER  OF  1E1S 

DELI  A- AMT  bE  TwFEttt  i  Vet.  S 

•-■■■•■-  =■  ■  ‘  ‘  •  "• 

RANCrx.RlSbtT--ETT£NTOF  A-PlANE  PLOT  IN  CYCLES  PER  POSITION  UNIT 

C  XAST  .VAST— COuRO  OF  Pf  FROM  «hKh  H*NGtX  *R*»i  r»fcV  MtASURtU 

READ  INPUT  TAPE  2  •  loo  1 .  f  r  1 A I  *4  1 A  i  •  »  U  I  **.*  1  *NJ 

1A0I  FORMAT  142k. 3F10.Sk 

i  . 

09  111  *«l»N 

'  --  - ===== 

==-•-- 

III  WIK1M* 

_  _  ' 

WRITE  0UfPUT  TAPE  3*3001*tYfAJ*X<*).«iM*«.*l*NJ 

3001  FORMAT  (UMuVAlOtS  OF  Y  .X  •  «/  i  3F2U*5  I  ) 

SNORM-g,‘j 

DO  100  Ul»* 

-V-  *  . 

P® 

100  5MORM-SNORI.4M1II 

jF|SMOftMH17.n2,lJS 

—  :--  -  - 
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GO  TO  110 
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11*  CONTINUE 
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$U*1«0*B 
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FM2-M2 
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== 

« 

«*rRST*»FMl-l.l»DY 

DO  20  J«1»M 
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»F(AR0I160.165.165 
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WRITE  TAPE  8»(vAlu«k.)  »L«l.igC* 

*1  WRITS  OUTPUT  TAPS  A* 1004* Y1 

1002  FORMAT  1AM  YfcwFM.St  1 

WRITE  OUTPUT  TAPE  3 . IUU 3 . t VALUI M »A-1 . lUO ) 

1003  FORMAT ( IH  .1UF12.5) 

_ PWI;1.5  _ _ 

60  50 

47  WRITE  OUTPUT  YAPS  1*1662. YK 

•RITE  Out'  Ut  TAPE  1.1003.  1 VALUIIU  .S-l *100 5 

Q 

r~^=  =f  jr 

PLOT-0.0 

- 

•  •  ij 

.  ..  . .  -GO  TO  5Q 

jgjpfi 

AB  PLOT-1,0 

.0  CfiWHlBJt 

II  CORTilHji 

REWIND  8 
! F i plot >  5 2* 
51  GO  TO  999 


51.52 


»  lFiWtTA-0.05)  55.5».54  ,-j 
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2000  FORMAT  1 1H  »22KQEt:A  TO  SMALL i DELTA*  .F10.3  )  _ _  _ _ 

GO  TO  999 

54  l«R  I  T£  OUTPUT  TAP£  3 .2001  .OELTA 

2001  FORMAT  (lHl»15HlNCKEMENTb  OF  .Flu. 3) 

WRITE  OUTPUT  I  APE  1.2002 

4002  FORMAT <  1M  ,10MLl«MfcAO«Y400»U**H78VOU3466709Qi4*4>6?8g8l2  345fc7 

_ llg0I234547ftg01234S67B9QI234547»9OI2H56789ul2S4KT89312l4567fc9u  l 

C  PAGE  MtAOER 
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DO  69  LINE  -1.100 _  _ 

ft£*D"TAPF  l*f VALUIS) 1001 
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_ _ 59  VSLUd  ) aAOSF I VALOIS 1 )  /SftQftH _ _  _ 

astep  »kfst 

DO  68  J*1 »100 

RNUM1  J>«VALUiJ>  _  _ _ 


IDIOT- 1 
04  $Tfcl*  *WTfP 

CLEY!  -  SWDELTA 
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l 
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■  Otffe  *  1.0 _ 

GO  TO  62 

66  IF  U  -1)  67.72.67 


_ 41-lgJLC.HlnG  i»yj_ 


71  PR1MT(J)»AlPhAi2AI 
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PRINT! j |  «  ALP.U INCHAR | 
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_ _ &Q  _ _ _ 

>2  SSTEP 
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604  CONTINUE 
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2003  format  (  ih  ,u,6h  »Iuuai*6h 
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69  CONTINUE 
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APPENDIX  B 


PROGRAMS -EVALUATION  OF  GAIN  IN  FREQUENCY  WAVE  NUMBER  SPACE 

The  programs  described  and  listed  in  this  appendix  were  used  to  generate  cross 
sections  of  patterns.  If  instruction  39  and  the  one  GO  TO  39  preceding  it  in  the  program 
3DFT  are  removed,  the  same  program  and  subroutines  can  generate  cross  sections  of 
gain  at  any  specified  frequencies  for  any  linear  filter  and  sum  processing  scheme. 

Form  of  Data  Output 

The  3DFT  program  and  associated  subroutines  compute  complex  gains  of  a 
processing  scheme  along  a  ray  in  wave  number  space  for  a  fixed  frequency.  The  gain 
is  evaluated  at  increments  of  DELTAK  in  wave  number  from  zero  to  NK*DELTAK.  For 
each  of  NFREQ  values  of  frequency,  FREQ(I)  =  1,. . .  ,NFREQ,  the  gain  is  computed 
along  NTHETA  rays  in  wave  number  space.  THETA(I),  I  =  1,. . .  ,NTHETA.  The  gain 
along  each  ray  at  each  frequency  can  be  output  in  several  forms.  The  choice  of  points 
to  evaluate  gains  and  the  choice  of  output  format  are  made  in  subroutines  UTIL  and/or 
OUTPU.  Thus,  these  can  be  modified  without  destroying  the  entire  program. 

In  the  current  form  of  OUTPU  several  options  are  available  once  the  complex 
gain  has  been  computed  along  a  ray  in  wave  number  at  a  fixed  frequency.  The  gain 
magnitude  and  phase  can  be  plotted  using  a  Stromberg  Carlson  4020  plotter.  The  log 
of  the  magnitude  of  gain  and  the  phase  can  be  plotted.  The  gain  or  log  of  the  gain  can 
be  printed  out  in  addition  to  being  plotted.  These  possible  outputs  corref  iond  to  a  con¬ 
trol  variable,  NCTRL  taking  the  values  1,  2,  3,  or  4,  respectively.  Two  variables, 
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ST01U  and  STOIL,  set  upper  and  lower  levels  at  which  gain  magnitude  or  log  of  the 
gain  magnitude  are  clipped  before  a  plot  is  made.  The  routines  KWKPLT,  STOID, 
CAMRAV,  and  PLTND  used  in  the  programs  are  routines  for  plotting  using  the 
Stromberg  Carlson  4020. 

Use  of  the  Program  (Fortran  IV) 

The  first  data  card  contains  NRUNS  in  Format  (15).  That  number  determines 
how  many  sets  of  data  are  to  be  processed.  Then  follow  NRUNS  sets  of  data.  Each  set 
is  composed  as  described  below. 

The  first  card  of  a  set  of  data  contains  NS,  NFP,  TFIRST,  DELTAT, 

(2I5,2F10. 5).  NS  is  the  number  of  seismometers  or  subarrays.  NFP  is  the  number 
of  time  samples  in  the  impulse  responses  of  the  filters  on  the  output  of  seismometers. 
TFIRST  is  the  time  before  which  all  impulse  responses  are  zero.  DELTAT  is  the  time 
between  adjacent  time  samples.  For  straight  sum  pattern  studies  one  can  use  NFP  =  1, 
TFIRST  =  0,  DELTAT  =  0. 

The  next  NS  cards  contain  Y(I),  X(I),  I  =  1,... , NS  in  the  format  (42X,2F10. 3). 

fH 

The  variables  Y(I),  X(I)  are  the  north  and  east  components  of  the  position  of  the  Iu 
point  in  the  array. 

In  general,  a  sequence  of  NS  groups  of  cards  must  follow.  The  I1-*1  group 
contains  W(I,J),  J  =  1,... ,  NFP  which  is  the  impulse  response,  starting  from  TFIRST 
of  the  I1*1  filter,  which  is  used  on  the  1^  channel.  The  format  is  (5F12.7).  For  pat¬ 
tern  studies  all  of  these  cards  a^e  omitted.  Instruction  39  and  the  GO  TO  39 
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instruction  set  all  filter  weights  to  zero  or  unity  as  required.  These  commands  must 
be  removed  to  do  other  than  straight  sum  pattern  studies. 

The  next  data  card  contains  NTHETA,  NFREQ,  NK,  DELTA K  in  format 
(3Ib,F10.  5).  These  variables  were  described  in  the  preceding  section. 

The  time  required  to  compute  and  plot  1000  points  along  a  ray  in  wave  number 
space  for  an  array  with  25  seismometers  has  been  less  than  .  5  minutes  using  an  IBM 


in^  !i£t  DIMENSIONAL  POURllN  TRANSFORM— -LACOSS  AU6UST  11.196S 

*  _  _ , _ 

DIMENSION  $T*Tul»24lI  j  »C5T&l  U4G1  i  »x  <60u  i « r  <60C  F*T  <SOi  »w</b»SO » 

C  NS  —  NUMbER  OF  Si  1  SMOMETEHS 

C  NFP— NUMBER  OF  POINTS  IN  TIME  FILTERS 

X  TR«iT^-Tf*e~cr M*si  output  from  filter  with  impulse  Amite 
j:  AT  TIM*  1680 

Z  WfcTAT— TJM|  BETWEEN  PTlTE*  SAMPLES  -  -  ^ - 

Z  xTTi  *¥<  l »— X“TNu  y  COOftO  OP  Ith  SllSMoMtTER  <am» 

C  TUI  — TIME  AT  P.«  SAMPLE  IN  TEMPORAL  IMPULSE  RESPONSES 

C  Wll»Jt  —  weight  OF  ITh  seismometer  AT  JTM  time  <1MPUL$E  RESPONSE) 

C - NTflT--<0«TR0L5  PROORAM  TERMINATION  1 VALUE  SET  «T  Wilt! 

^  C  TRIVIAL  CXANWtS  OF  DIMENSION  STATEMENTS  WOULD  ALLOW  MORE  ON  LESS 
_  SEISMOMETERS  AND/OR  PlLTtR  POINTS  IF  NtMONT  Is  AvAIU#LE 

~"C - TiXIlffl^EAL  “PiRT'OF  TtsPofisr  Transform 

C  GAIN 1--1M AO  INARY  PART  OF  RESPONSE  TRANSFORM 

C  NRUNS--NUMBER  OF  DATA  DEC AS  SUttM I T T F D  _ 

MHITC"14t7®T 

tWI  THHFF  DIMENSIONAL  FOUR  I  IN  THAMSFORM-- LACOS5  AUG31*196S| 


iw  SN<UL(  il  •NlNUHLlAI 

7M=:=m*'U  f»im»gAOt  , _ _ _ 

T?  SNT BliTM-Sfi f  dL  I  I  *■1/00  ) 

DO  14  i-l.UOl 

14  CSTbH  I  PS-NTBH  I  »6CiUI  _____ _ _ 

L  “  D^UT-iiOifiMil 

t  *SN  TAUI-ISOO  j 

C  t  BIN  COS  TABLE  IS  COMPLETE  .  _  - -  - 

C  INITIALIZE  DATA 

RE AD( S  «  S99  > NRUNS _  _ 

-5WF<8U0Tn  5 1  "  .  - 

i»  5U0  kfcfck>}»NNU*$ 

_=l==  READ 3 StSOff I  _NS*NFPt TF 1  «S T» OELTAT _ -  _ _ _ 

600  FORMA?  <21**2F10*S) 

WRITE  <6*7011  NS*NFP*TF1RSTiDELTAT 

701  FORMAT  |4HQNSyl  3*  I  OX  *  hHNF  P  •  |  3jJ  L»X»7HTF1RST-F9,S*1UX  , 7HDEL TAT»F9. 

READ  I5.4C11  UU  I  *X<  1 1 .  l»l*NS» 

FORMAT  f«X«2F10*S  I 

_  (YUlfAUUl^tjRil  ;  ~ 

702  FORMAT  <  36HUSE 1 SMOME  TtR  LOCATIONS  N-S.E-W  IN  XM(2F20*3»10X*2F2U« 
DO  30  1  •  1  *NFP 

_ FLIjlI  _ _ _ 

30  fm*TPiRST40CLTAT*lpLl-l,i 
00  40  Witt 

"read<5»602T  (W< 1 *J| *J«1.NFP| 

60?  FORMAT ( 5F1 2*  7  I 


m«4NTITE  1 6 *?0A1  Ft  l*J| I  *Ul  .0*1  •  NFP »_ 
PORMAT(20*4OlMPDLSE  RESPONSE  UFfti 

~  J*L*GZ Jl_ _ _ _ _ - 

C  ENO  OF  INITIALIZE 


re/IIOFlOrSTI 


jlgllpgF^W  WBDtte  1 5<Wf*NPU$ 

CALL  UTlUXAtTKtfREOtfiAlMRtDAiNJtNFLAStil 
GA1NR-C. 

GAIN  1 *0* 

_  DO  200  I ■ 1 » NFP 

m  »  tot  o-i. ns 

L  :— -  AR€^RR*XlJl*YR*VT4I+P«C0*Tfll 

iiRfr*AMDP|AR6»i«l  .  _ 1_ 

|F  < ARG 160*66  *66 
60  ARG»ARG*1 • 

6S  ARG  •24O0*«ARG 

=  2  _  K* I F  1*1  aROT*! 

BBE.  '**TIO«AMufr|ARO»i«  1 _ _ _ 

COSVAL •CSTdL I M  *RA  T 10#  <  CST  tiL ( XX I -LS  TdL ( A  I » 
S|NVAL»SNTBLIM*RAT  lUMSNTbLtkXI-SNTbLIX)  I 

_ oA1NR»GA1NR*W<  J.  WlUSVAl.  _ _ 

i=-  _JY60aAINI»flA|NL+WI  J*  1  t*SlNVAL 

WIfeWt«£FfREOrOAlMfc,OAiNi,NfLA«,li 

_  __  ftQ  TO  NS  . _ --  ■■ _ _ 

S00  CONTINUE 
CALL  PLTND 
-  STOP 
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v?ti.t)w«T*«f *64i«H**6Ai4i«wi4t,jNWTj _ 

$£PT  2. 1<»6S" 

ROUTINE  SUPPLIES  XI.(HtF  TO  FUUNltR  T RANiF okME*  AND 
O’JlPUTS  VALUES  OF  TRANSFORMS 


D|W£NS’0«  TH£T4<20I.FR£OUUI«STOiUWOOI.St02aOOO» 
m  TO  (1*1  WOT  *  T*WUT 
.  1  60  10  l> * 50*50 1*NFLAG 


INITIALIZE  UN  FIRST  ENTRY 
5  WHITE  (6.7UOJ 


TOO  FORMA  T  (33HUUSEU  Ulkllr  RKUGKAM  SENT  2.1V6S) 

MAS  15.400 1  NTMETA.NFftEu.N*.. DELTA*. 
m  400  fORMATiJIS.f 10*51 

Z, _ TWITE  14*7011  NIHETA.NFRFQ.  Hk,  OELTAx _ 

701  FORMAT  ( 1  OX  .  7hNThET  A«1  5  .  YX  *t>MNF  RtU«  I  5  *  YX*  inNx«  »&»*X*7riDELIAx«HU»5l 
C  HTHETA — NUMdER  OF  VALUES  OF  CLOCRWlSt  ANyLtS  Fwu*  NORTH  IImEIAI 

_ C _ jHA_T_ARt  TO  Ot  UStO _ _ 

c  nfreo-- mmt*  of  freu^nc*  if  or  freoi  to  «  usto 

fe=slc.^-.  IB  —  NUMBER  OF  RAO J  At  VALUES  OF  WAVE  NUMBER  10  BE  USEO 

C  DELTAS— DIFFERENCE  dtTWEEH  ADJACENT  PAQIAl  VALUES  Of  W*VE  NyMBEf* 

READ  (5*60!)  (T-<ETAt  i  |  »I>l,NTHkTA) 

READ  (5»6U1>  <  FREOt 1 ) il'ltNFREu) 

401  F0RMAT{0F1U,>) _ 

tel"-  ■_  WRITE  (4*702  1  (THETA!  1) 1 1-1  tMTKETA I 
^  702  FORMAT i 1 6H0 VALUE S  OF  TrtETA/l lCFiu.S) ) 

K,  ^  . WRITE  14t?V3UFRfc«fiUl*l*MFRE0i 

703  FORMAT ( 1  SHU VALUES  OF  FRkU/ i  10F 10.4 ) ) 

ITHETA.1 


_ _ iEBEfisl _ 

IR«1 

WFLA6«l 

n=--^.  --  h^eAs*? _  ...  _ _ 

C  INITIALISATION  IS  COMPLETE 


C 


S5.55.60 

—  %«c  TEMP^THElM  (  I  THETA  » 

CTHETA.TOSf  TEMPT 


STHkTA»S!N( TEMP) 


60  STZEX*  (FLOAT  I  I  M~1  •  I  *DELT  A*, 

#K«$THEIA*$IZEK 
T*«CTHETA«5UEK 
F»FREQ( IFREO) 

ffTTOffTT- 


C  XX.YX.F  HAVE  BEEN  GENERATED 

C 

LT  C  'Eimv^POIhT  FOh  IS0UT42 
100  STOH  UMGAINN 
ST02(Jxi»OAf*I 
fTTNX*U  i  llu.liv.lTu 


110  CALL  0UT°UT (STUl .STU2.HFLA0.NX.0ELTAM 


MF  LAG*  2 

~7t ( WTRITj^TThIT a!~II5 . i 15 . i*o 

115  IMN*«t«-lFRfcUl  120*120.130 

;.  _  UO  mflaoo _ 

RTTuRn 


130  IX«1 

1THETA.1 

-  TFREO-IFRECU 
RETURN 
140  U«I 

TThETa^TThTTa,  :  ~ 
RETURN 
ISO  I  X *  I  X ♦ 1 
>4V/  RETURN 

^  tWL  — - 
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SUbWOUl  i**t  out Pul  i  aTvil  *3  F02  **FlAu»Ni  »UtL  f  *x.  I 
C  $tPT  iil'tb* 

Itu  Hft’fbb 

L  PLOlb  TRANSFORMS  OF  3f»IM  FUNCTIONS 

dj«e»$ion  s*eiftfcOv*wT«iUvse> 

n IMF NS I  ON  XCOOR ) 

GO  io  I  1  ■%!  tMFi.AC> 
l  BfeAP  I5.6JUI  hvlRi.  tbTQiv#:  Tvit 
too.  FORMAT  e  {  t,2F  !U«bJ 

«Sll£  *6.7iCj  HitEt  »5T0iy  *51014. 

700  f  ORMAl  !  lSHvOUfPul  CONT  Rui./ 7n  N(.  «  *«L  *  i  N  •  R*  #Nrr->  v  01  V=P  ,  «  x  *  &.\n  |  v  1 

IFlO.Si 

C*LU  S 1 0 1  Lf  V  l  £  '-•  Mu  A  i  H  Afto  PHAGE  rUi’i.i 
C*LL 

SO  IVcv  P«jtNF 

iO>0  *tO*OR4U» J-l » «0£l! Ak 
C  1KT4L-* «IFOI  tuNIRUL 

t  I— PLOT  GA1S,P«A5fc 

C  2— plot  LOOuAIN.PmASE 

C  3—1  «N0  PRINT  5A**E 

C  «*--£  A  NO  PRINT  SAVL 

C  StOlU— UPPcR  CLIP  LEVEL  FO:x  bTul 

C  5T0U--i.C*€R  CLiP  LEVEL  FOP  ST01 

C_  max.njn  Fqe  $TC2  ARE  »P I  AND  -PJ 
5  DO  iU  |».  t Hr. 


_i£MPiiiaiiii**z*aiu2Ui«*2 

GAIN.SORUT'Hki 

TEWP^S?v<(i)  .  .  _  _ 

PHA5E**TAN^|TtMF^,U*Mi, 

$7011 IT-OAJN 

IFll.FO.l  »  wC«ST0H1J _ 

STOllI  »*$T01in/0C 
Ifr  $T02ft7»PHA$L 
-  .  _  VRlfElfr  «  702 »_UC_ 

702  FOR*aU23mv/2Ehu  WAVE ‘nuMHEK  GA|N«F2u,7) 

GO  TO  t*U#20#3u*2U1*NCIHL 
20  DO  25  1*1. NA 
25  T«MP*5r&mi 

5T0H  1  }*2tf*«ALOG4  TEMP) 

_  6C  TO  <*S»*S»10»3ol»NCIHL  a _ 

30  WR  (t»»7vj>  I  5  T  0 1  (  tl  »  3 1  Oi  !  I  J  » l 3 1 • NA  } 

701  FORMAT  (///42F’,y«5t9*,2‘  iu.b  lyA«^lUt3i  VA  .2H0.5  » 1 


AO  CONTINUE 

00~ 50  !*!#«* 

**4.$T01tl)«*C«»fOlU»  ST01  ( I  |a&f010 
~X-  SO  lfi5T0|finteS|glkj  5T01U)»5I01L  _ 

STvlU  j*5T01L 

CALL  AWAFtl IXCOUH.SlOl  .NA.1BM  «AVE  NuMoLl*  .IBM  GAIN 

_ 1  _ I 


£=rirr  . 

GALL  AIAPlT tALUOK.bltU  *NA«2IH  RAVE  NWNoLR 

•  liM 

phase 

S5£g|s?-*’- 

l  f 

HFLAW 

RETURN 


END 
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A  PPEND1X  C 


USEFUL  PATTERN  SYMMETRIES  AND  INVARIANCES 

A  comparison  between  members  of  the  same  family  of  geometries  or  between 
different  families  is  made  a  bit  less  difficult  when  certain  symmetries  of  patterns  are 
noted  and  invariances  with  certain  changes  of  geometry  are  taken  into  account. 

For  any  geometry  whatsoever,  it  is  true  that 


i H(k  ,k  )|  = 

i  x  y 


!»/. 

Ill\ 


-k  )\ 

Y  ' 


This  follows  from  substitution  into  (1)  and  noting  that  a  number  and  its  complex  conju¬ 
gate  have  the  same  magnitude. 

Suppose  that  a  geometry  is  transformed  into  itself  by  a  rotation  of  a  radians. 
Then  the  pattern  of  that  array  is  transformed  into  itself  by  a  rotation  of  -n  radians. 
This  is  best  seen  by  using  polar  coordinates  k,  3  for  wave  number  and  r^,  jr  for  the 
location  of  the  q  n  element  of  the  array.  Using  this  notation,  the  expression  for  the 
pattern  becomes 

1  - 

II(h,9)  =  ^  (j  exp  [ -i2~kr  cos  (9-t  )1  • 
q=l 


Reproduction  of  the  array  by  rotation  of  a  radians  means  that 


i  Ni 

rl(k.r)  =  —  exp  {  -i  2”kr  cos  [ ~  —  (c  +c)]  1 
N  “  q  q 

q=  l  ' 


But 


cost?  -ft  +a)]  =  cos[(3-a)-'Ll  . 

4  4 

Using  this  fact,  one  obtains 

H(k,  9)  =  H(k,8-a) 

as  as  asserted. 

The  facts  derived  above  imply  the  following  results.  Let  a  be  the  smallest 
rotation  which  transforms  an  array  into  itself.  If  rr  is  an  integral  multiple  of  a,  then 
the  pattern  for  that  array  will  repeat  itself  every  a  radians.  II  n  is  not  an  integral 
multiple  of  a,  then  the  pattern  for  the  array  will  repeat  itself  every  a/2  radians.  The 
first  of  these  results  is  true  directly  from  the  relation  H(k,  9)  =  H(k,  8-a)  derived 

,-ky)|  for  any 

geometry  at  all. 

Finally,  consider  a  geometry  in  which  there  is  a  ring  of  elements.  Suppose 
that  elements  are  uniformly  spaced  on  that  ring  and  that  one  element  is  found  due 
north  of  the  array  center.  A  rotation  of  that  ring  by  rr  radians  leaves  the  pattern  un¬ 
changed  in  the  east  and  west  directions  in  wave  number  space. 

This  fact  is  seen  as  follows.  The  function  H(k,8)  is  a  one -dimensional  transform 
of  the  array  geometry  projected  onto  a  line  in  the  9  direction.  The  rotation  of  a  single 
ring  as  described  above  leaves  the  projection  of  the  array  onto  a  line  in  the  east -west 
direction  unchanged.  Thu®,  H(k,  9)  is  not  modified  if  9  is  the  angle  which  idenl  fies 


above.  The  second  requires  use  of  the  fact  that  |  H(kx,k^)  |  =  |H(-k^ 
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either  east  or  west.  As  a  direct  result  of  this  invariance,  one  can  conclude  that 
patterns  for  geometries  such  as  are  shown  in  Fig.  6  must  be  the  same  in  the  east  or 
west  directions  in  wave  number  space.  In  the  particular  case,  the  same  is  true  in 
directions  an  integral  number  of  tt/5  radians  from  east.  This  invariance  with  rotation 
can  be  useful  when  trying  to  compare  geometries  with  an  odd  number  of  elements  per 
ring  with  and  without  rotation  as  shown  in  Fig.  6.  This  invariance  is  not  relevant 
when  consideriqg  geometries  having  an  even  number  of  elements  per  ring  such  as 
those  shown  in  Figure  7. 
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APPENDIX  D 


PROGRAM-SENSITIVITY  SUBROUTINE 

Plots  of  sensitivity  of  patterns  with  radii  of  rings  have  been  obtained  using  the 
subroutine  listed  in  this  appendix.  The  routine  is  called  by  3DFT  which  is  listed  in 
Appendix  B.  The  pattern  along  a  ray  in  wave  number  and  the  sensitivity  along  that  ray 
for  each  ring  are  computed.  This  information  is  plotted  out.  This  can  be  done  for 
several  directions  in  wave  number  space. 

Use  of  the  Sensitivity  Programs 

3DFT  must  be  submitted  as  the  main  program  to  be  used  with  the  sensitivity 
routine .  The  output  from  the  program  as  written  requires  the  Stromberg  Carlson  4020 
plotter.  If  this  is  not  available,  the  basic  program  remains  valid  but  output  instructions 
must  be  changed. 

The  first  data  card  contains  NRUNS  in  format  (15).  This  is  the  same  variable 
described  in  Appendix  B  and  indicates  the  numbe?  of  sets  of  data  to  be  processed. 

Each  set  of  data  is  constituted  as  follows.  The  first  card  contains  NS,  NFP, 
TFIRST,  DELTAT  in  formac  (2I5.2F10. 5)  with  NFP  =  1,  TFIRST  =  0,  DELTAT  =  0, 
and  NS  the  numbe-  of  seismometers  or  subarrays.  The  next  NS  cards  contain  north 
and  east  coordinates  of  subarrays  or  seismometers  in  a  format  fc2X,2F  1(1. 3).  The  next 
data  card  contains  NTHETA,  NK,  NR,  NPR,  DELTA K  in  format  (415,  F 10. 5),  NTHETA, 
NIC,  aud  DELTAK  have  the  same  meaning  as  in  Appendix  B.  It  is  assumed  that  the 
array  has  at  least  NR  rings  of  NPR  evenly  spaced  seismometers  or  subarrays.  The 
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sensitivity  is  to  be  computed  for  the  radii  of  these  NR  rings.  The  next  group  of  data  is 
RADIUS(I),  I  =  1,. . .  ,NR  in  format  (8F10.5).  These  are  radii  of  rings  of  interest. 
Finally,  ANGLE(I),  I  =  1, . . .  ,NR  in  format  (8F10. 5)  are  read.  ANGLE(I)  is  the 
orientation  in  radians  of  one  element  on  the  I1*1  ring,  which  has  radius  RADIUS(I). 
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SU0ROUT  1 NE  UT1UXH.YI  .  F  *GA I  NR  »GA 1 N I » NF L AG »  INuuT  > 

C  OCT  30*1965 

C  ROUTINE  SUPPLIES  UiTK.F  TO  FQURIEn  TRANSFORMER*  CUMPUIES 

~c  ^  smitmrfi?  mm&&m  *tmnv#L  *m  «aou  of  rings  of 

^  ARRAY  •  PLOT  A  LEASH  IV)  fY  AMO  TRAM  FORM  MAftMl  Ifi$€ 

C  A  RAY  \k  *AVg->U*«R  $PA« _ _ , _ '  '  ~  ~  — 

"BTl^fiNS |ON  ThEIA(2Q>  .ST01UOQQ)  »GTQ2l  10UU>*S£NS110UO>  *RAC(U5I10> 
DIMENSION  ANGLE < 10) • PS t ( 23) *XCOOR! 1000) 

GO  KM  1  *10U)  *INUUT 

Hi  TTW  T0lf«50.501  »Wl*4  '  '  ■ '  ' 

OMFll&TfcMIRY 

_  *  WUJEU^TO*  _  _  ■■  ■  ■  _  ■-_  .-  -  -  _ 

"fB6“?0RMAft33H?USlB  uTFlITy  PROGRAM  UCT  30*1965) 

READ!  5  .600) N THETA ,Nl. NR .NPR.DEL  TAi 
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ADJACENT  RINGS  ROTATED  RELATIVE  TO  EACH  OTHER 


NO  ROTATION  BETWEEN  ADJACENT  RINGS 

Figure  5  Array  geometries  having  seven  subarrays  per  ring 
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NO  ROTATION  BETWEEN  ADJACENT  RINGS 


Figure  6  Array  geometries  having  five  subarrays  per  ring. 
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(A)  Array  geometry. 


Figure  9  Geometry  and  patterns  for  an  array  having  seven  subarrays  per 
concentric  ring  and  geometric  spacing  of  rings. 
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Pattern  (contour  map) 
Figure  9  Continued. 
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(A)  Array  geometry. 


Figure  10 


Geometry  and  patterns  for  an  array  having  seven  subarrays  per 
concentric  ring  and  geometric  spacing  "f  rings. 
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Figure  11 
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(A)  Array  geometry. 


Geometry  and  patterns  for  an  array  having  five  subarrays  per 
concentric  ring  and  geometric  spacings  of  rings. 
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Figure  11  Continue 
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Figure  12 


(A)  Array  geometry. 


Geometry  and  patterns  for  an  array  having  five  subarrays  per 
concentric  ring  and  arithmetic  spacing  of  rings. 
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Figure  12  Continued. 
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(C)  Pattern  (cross  sections). 
Figure  12  Continued. 
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(B)  Pattern  (cross  sections). 


Figure  14  Geometry  and  patterns  for  an  array  having  four  subarrays  per 
concentric  ring. 
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Figure  17 


Patterns  (cross  sections)  for  both  arrays  shown  in  Figure  15 
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Figure  20  Patterns  (cross  sections  in  the  direction  18°  north  of  east  in  wave 
number  space)  showing  some  effects  of  reducing  the  25  km  ring  of 
Figure  12A  to  20  km . 
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Pattern  (cross  section)  and  sensitivity  functions  for  an  array  geometry 
which  would  be  obtained  from  that  of  Figure  11A  by  increasing  the 
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Figure  23 


Pattern  (cross  sections  in  the  direction  9°  north  of  east  in  wave  number 
space)  showing  some  effects  of  increasing  the  55  km  ring  of  Figure  11A 
to  57  km . 
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Figure  24 
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Pattern  (cross  section)  and  sensitivity  functions  for  the  array  geometry 
shown  in  Figure  9A . 
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25  Pattern  (cross  section)  and  sensitivity  functions  for  the  array  geometry 
shown  in  Figure  10A  . 
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